Electromagnetic Tunneling through a Single-Negative Slab 
Paired with a Double-Positive Bi-Layer 
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We show that resonant tunneling of electromagnetic fields can occur through a three-layer struc- 
ture composed of a single-negative (i.e., either negative-permittivity or negative-permeability) slab 
paired with bi- layer made of double-positive (i.e., positive permittivity and permeability) media. In 
particular, one of the two double-positive media can be chosen arbitrarily (even vacuum), while the 
other may exhibit extreme (either near-zero or very high) permittivity/permeability values. Our re- 
sults on this counterintuitive tunneling phenomenon also demonstrate the possibility of synthesizing 
double-positive slabs that effectively exhibit single-negative-like wave-impedance properties within 
a moderately wide frequency range. 

PACS numbers: 42.25.Bs, 41.20.Jb, 78.20.Ci 



The past decade has witnessed a growing interest in 
the study and applications of extraordinary transmis- 
sion of electromagnetic (EM) fields through opaque me- 
dia, such as metals or, more in general, single-negative 
(SNG) materials characterized by negative permittivity 
or permeability, and hence evanescent fields. With ref- 
erence to metallic films, surface-plasmon-induced trans- 
parency was demonstrated theoretically and experimen- 
tally in Ref. fl and a large body of results are available 
in connection with extraordinary transmission through 
subwavelength apertures (see, e.g., Ref. d for a recent 
review). In connection with general SNG media, it was 
shown in Ref. 3 that, in spite of their inherent opacity, 
slabs of (homogeneous, isotropic) epsilon-negative (ENG) 
and mu-negative (MNG) media can give rise to resonant 
tunneling phenomena (with total transmission and zero 
phase-delay) when paired as a bi-layer. Starting from 
this basic configuration, several extensions and general- 
izations have been investigated (see, e.g., Refs. [j-Qj]). 

Moreover, multi-layer metallo-dielectric structures 
have been extensively studied in connection with near- 
field sub- wavelength imaging (see, e.g., Refs. [r3 - [l5h . 
generating a growing interest in the study of resonant 
tunneling of evanescent waves in configurations featur- 
ing SNG media paired with double-positive (i.e., positive 
permittivty and permeability - DPS) materials. Inter- 
estingly, while it is impossible to achieve perfect trans- 
parency by pairing a single SNG slab with an arbitrary 
DPS one, anomalous tunneling phenomena have been 
demonstrated in configurations featuring an SNG slab 
sandwiched between suitable impedance-matching DPS 
layersi 16 ' 17 In particular, the results in Ref. per- 
taining to a symmetrical DPS-ENG-DPS three-layer and 
experimentally validated at microwave frequencies, have 
been extended in Ref. 18 to more general (e.g., opti- 
cal, quantum-mechanical) tunnel barriers. Against the 
above background, we propose here another counterintu- 
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FIG. 1. (Color online) Problem geometry (details in the text). 



itive resonant tunneling mechanism which entails pairing 
an SNG slab with a bi-layer made of DPS media. Unlike 
the configuration in Ref. [la, where the SNG slab is sym- 
metrically sandwiched between two identical DPS layers, 
in our case the DPS layers are different and paired at one 
side only. This also allows direct comparisons with con- 
figurations involving paired ENG-MNG media, with the 
DPS bi-layer playing the role of an effective SNG slab. 

The problem geometry is illustrated in the Cartesian 
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(x, y, z) coordinate system of Fig. Q] Without loss of gen- 
erality, we consider a homogeneous, isotropic ENG slab 
with relative permittivity E\ < and thickness d\ paired 
with a bi-layer composed of homogeneous, isotropic DPS 
media with parameters £2 > 0, di and £3 > 0, g?3, respec- 
tively, all embedded in vacuum (e = 1) and assumed 
as non-magnetic (i.e., /i = 1), under time-harmonic 
[exp(— iujt)] plane-wave illumination. 

Assuming normal incidence, for the ideal case of loss- 
less media, we can derive via straightforward transfer- 
matrix algebra^ the general resonance condition for total 
transmission, by zeroing the reflection coefficient, viz., 



i (1 - £1) v^2£3 T i + * (1 - £2) s/-£\£zt<i 
+ i (1 - £3) v / -£i£ 2 t 3 + v / £i(£ 2 - ei)nr 2 
+ V^2 (£3 - ei)riT3 + V-£i (£3 - £2)t 2 t 3 
+ i (£i£ 3 - £ 2 ) riT2T 3 = 0, 



(1) 



where t\ = tanh (k-^— e\di) and r 2j 3 = tan (^^£2^3(^2,3) , 
with k — lu/c = 2tt/\ denoting the vacuum wavenum- 
ber, and c and A the corresponding wavespeed and wave- 
length. Zeroing the real part of ([lj, we obtain 
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y/£3 (£2 - £i)n + s/-E\ (e 3 - £ 2 ) r 3 
which, enforced in the imaginary part of ([1]), yields 



(2) 



T3 = ±1 



£3 (1 - £1) (£2 - £1) rf 



y £1 (1 - £3) (£3 - £2) + (£3 - £1) (£3£i - £2) t\ ' 

(3) 

Noting that the numerator of the square-root argument 
in ([3]) is always positive (since £1 < and £ 2 ,3 > 0), real 
solutions exist iff 

£1 (1 - £3) (£3 - £2) + (£3 - £1) (£s£i - £2) t\ > 0. (4) 

The above inequality is quadratic in £3 (with positive 
discriminant), and yields the conditions 



£3 < £3a, or £3 > £ 3b , 



(5) 



with £3 a and £36 denoting the two (positive) roots of 
the denominator of the square- root argument in ([3]). 
It can be readily shown from (QJ that £30 < £2 and 
£36 = £2/£3a > £2, from which it emerges the impossi- 
bility of achieving total transmission by pairing an ENG 
slab with a single DPS slab in vacuum, consistent with 
Ref. [HI. Conversely, in the presence of a DPS bi-layer, 
for a given ENG slab with parameters £1, d\, it is always 
possible [from ©-([S])] to find suitable parameters for the 
DPS layers so as to fulfill the total-transmission condi- 
tion in (JlJ at a given frequency. More specifically, the 
parameter e 2 turns out to be completely arbitrary (but 
positive), while £3, e2 3 , and d 2 are derived proceeding 
backwards from (O, (|3|), and ([2]), respectively. 

A few general considerations are in order. First, at 
variance with the symmetrical DPS-ENG-DPS configu- 
ration in Ref. 16, no critical thickness value exists for 



our asymmetrical ENG-DPS-DPS configuration in order 
to exhibit EM tunneling. In fact, there always exist four 
distinct classes of solutions for total transmission, cor- 
responding to the possible combinations of the choices 
in ([3]) and ([5]). We found that, in the long- wavelength 
limit d n /Xo — > 0, n = 1,2,3, all these solutions tend 
towards the effective-medium-theory prediction obtained 
from (TT]) by neglecting the multiple-scattering terms (i.e., 
those containing multiple r-terms), 



e%dx + e 2 d 2 + £3^3 = d\ 



(6) 



It can be shown that, for an increasing opacity of the 
ENG medium (|ei| >• 1), the end values in ([3]) tend 
to exhibit extreme values, so that the outermost DPS 
slab must have either near-zero {e^ a <C 1) or very-high 
(£36 3* 1) permittivity. From ([3]), it is readily verified 
that the solutions associated to the end values £3 = £3 a or 
£3 = £3b in © correspond to a quarter-wavelength (plus 
half-wavelength periodicities) size for the outermost DPS 
slab, which (recalling its wave-impedance-transformation 
properties-^) is also representative of a different config- 
uration featuring an ENG-DPS bi-layer terminated with 
a DPS half-space (i.e., d% — > 00). Moreover, we highlight 
that the complete arbitrariness in the choice of £2 rep- 
resents an important degree of freedom in the proposed 
configuration. In particular, choosing £2 = 1, we obtain 
another interesting tunneling configuration featuring an 
ENG slab paired with a DPS slab via a separating vac- 
uum layer. Finally, it can be shown that enforcing a 
purely real (i.e., ±1) transmission coefficient results in 



H - £ i £ 2 - £i£3 + £2£3j = (1 - £1) (£2 - £1) 
'£1 
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r(£3-£2) (£ 3 -l) + (£l-£ 3 ) (£l£3-£2) 



(7) 



i.e., a quadratic equation in £3, whose solutions, enforced 
in (j4]), yield the constraint 



e?(l-ei)(ei -s 2 f^ > 0, 



(8) 



which is clearly impossible to fulfill with S\ < and £2 > 
0, thereby implying that no tunneling with zero phase- 
delay is possible in our proposed configuration. 

As a first example, we consider an ideal lossless con- 
figuration featuring an ENG slab with £1 = —3 and 
d\ = 0.1 Ao (here and henceforth, the subscript iden- 
tifies resonant frequency /wavelength quantities). We ar- 
bitrarily select £2 = 2.5, and readily derive from ©-([S]) 
the remaining parameters: £3 = £ 3 b = 16, c?3 = 0.0624Ao 
[choosing the positive determination for T3 in ©], and 
d'2 = 0.236Ao. For a normally-incident (along the posi- 
tive ^-direction) plane- wave excitation, Fig. [5] illustrates 
the resonant electric and magnetic field (normalized) lon- 
gitudinal distributions, from which it can be observed the 
total-transmission effect, achieved via a growing evanes- 
cent wave in the ENG layer and a standing wave in the 
DPS bi-layer. The electric field peaks towards the center 
of the middle DPS layer, and the magnetic field peaks 
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FIG. 2. (Color online) Resonant electric (solid) and magnetic 
(dashed) field intensity distributions (normalized by the inci- 
dent fields) for a configuration as in Fig. [1] with ej = —3, 
di = O.lAo, e 2 = 2.5, d 2 = 0.236A , £ 3 = 16, d 3 = 0.0624A . 
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FIG. 4. (Color online) Black-solid: Transmittance as a func- 
tion of the frequency for the parameter configuration in Fig. 
[2J but using for the ENG slab the dispersive model in (|9]) 
with cjpi = 2wo, 71 = 3.75 • W~ 3 ui p i (i.e., Re[e3(aio)] ~ —3), 
and a loss-tangent of 10 -3 for the DPS layers. Also shown 
are the responses obtained for an increased opacity of the 
ENG slab (u pl = 10.05u and 71 = 9.8 • 10~ 4 lj p1 in ©, i.e., 
Re[ei(cj )] « -100, and di = 0.01A ), for e 2 = 12 (with loss- 
tangent=10 -3 ) and d 2 = 0.119Ao (red-dashed), and £2 = 1 
and d 2 = 0.472Ao (blue-dotted), with the high-permittivity 
outermost DPS layer described by the model in QQ) , with 
£300 = 4, A 3 = 3.081cj 3 , and u) 3 = 1.091w , 73 = 9.26- 10" 4 cj 3 
(i.e., Re[£3(tJo)] = 63.7), d 3 = 0.036Ao, and ujz — 1.13a;o, 
73 = 1.33 • lQ- 3 u) 3 (i.e., Re[£ 3 (^o)] = 47.8), d 3 = 0.0313A , 
respectively. The inset illustrates the influence of losses in 
the ENG material on the peak transmittance (as a function 
of the loss-tangent at resonance). 
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FIG. 3. (Color online) Reflection coefficient magnitude (a) 
and phase (b) of the isolated DPS bi-layer in Fig. [2] (blue- 
solid) as a function of the frequency, compared with the re- 
sponse (red-dashed) of an effective MNG slab with e e = 1, 
= —0.33. 



at its boundaries, similar to a Fabry-Perot etalon, al- 
though we consider here the presence of an inherently 
opaque material. It is insightful to compare the above 
phenomenon with the tunneling occurring in an ENG- 
MNG pair. From the ENG-MNG matching conditions 
[cf. Eq. (8) in Ref. yf, we can straightforwardly de- 
rive the constitutive parameters e e = 1 and fi e = —0.33 
of the required MNG slab of same thickness (d^ + cZg) 
as the above DPS bi-layer to "compensate" the given 
ENG slab. Figure [3] compares the reflection-coefficient 
responses of the isolated DPS bi-layer and such MNG slab 
(free-standing in vacuum) as a function of frequency. As 
expected, both responses result in large reflectivity (due 
to the opacity of an MNG stand-alone layer) and perfectly 



match at the resonance frequency. In particular, the de- 
signed DPS bi-layer is capable of providing the capaci- 
tive input impedance required to resonate with the ENG 
slab, ensuring total transmission. The reflection coeffi- 
cients in Fig. |3] agree reasonably well within a moderate 
bandwidth (maximum difference of ~ 0.5% in magnitude 
and ~ 7r/12 in phase, over a 10% bandwidth), ensuring 
that the designed DPS bi-layer may effectively replace an 
MNG layer for a variety of applications. 

Next, we move on to assessing the frequency and angu- 
lar dependence, as well as the sensitivity to polarization 
and losses, of the above tunneling phenomenon. Figures 
3]and[n]show the transmittance as a function of frequency 
and incidence angle (from the a;— axis, for both P and S 
polarizations), respectively, for three representative pa- 
rameter configurations. We start from the same con- 
figuration in Fig. [21 but now assuming a more realistic 
Drude-type dispersive, lossy model for the ENG medium, 



e x (w) = 1 - 



lo (ui + Z7l) ' 



(9) 



with the plasma angular frequency u p i and the damp- 
ing coefficient 71 chosen so as to ensure Re[ei(o;o)] ~ 
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FIG. 5. (Color online) As in Fig. [4] but as a function of the 
incidence angle at the resonant frequency, for the P (a) and 
S (b) polarizations. 



—3 (with a loss-tangent ~ 10 -2 ), and a non-dispersive 
slightly-lossy (loss-tangent=10~ 3 ) model for the DPS 
layers. From Fig. @] (black-solid curve), a rather 
broad resonance is observed in this case, with a peak- 
transmittance of nearly 97%, and a slow decay for higher 
frequencies attributable to the increasing transparency 
(approaching the plasma frequency) of the ENG slab. 
Also the angular response (see Fig. [5]) turns out not to 
be very selective, especially for P-polarization. As previ- 
ously highlighted, tunneling effects may be obtained, in 
principle, for arbitrary choices of the ENG slab parame- 
ters and the permittivity of the middle DPS layer. For 
an increased opacity of the ENG medium, as anticipated, 
the constitutive parameters of the outermost DPS layer 
tend to exhibit extreme values. For instance, consider- 
ing an ENG slab with e\ = —100 and d\ = O.OlAn, and 
£2 = 12, we obtain [choosing the larger end- value £3;, in 
©] £3 = 63.7. In Ref. [TH, such (positive/negative) high- 
permittivity media were successfully synthesized at mi- 
crowave frequencies with resonant inclusions. Assuming 
for the outermost DPS medium a Lorentz-type disper- 
sive, lossy model, 

A 2 

£3 M = £3oo - -3 27T . ( 10 ) 

with the parameters (given in Fig. |4]caption) tuned so as 
to ensure the required real-part at the given frequency, 



we observe in Fig. [4] (red-dashed curve) a narrower band- 
width (with peak-transmittance of nearly 85% due to the 
larger sensitivity to losses, and a smaller transmission 
peak attributable to dispersion effects), and from Fig. [5] 
a flatter angular response (especially for P-polarization) 
which is a direct consequence of the increased permittiv- 
ity values.— The resonant field distributions, not shown 
here for brevity, are qualitatively similar to those in 
Fig. [2] with intensity enhancements of nearly a factor 
60. Also shown in Figs. 0] and [5] (blue-dotted curves) 
are the responses obtained for the same ENG slab, but 
choosing £2 = 1, i.e., an ENG-vacuum-DPS configura- 
tion. While the frequency response and the field distri- 
butions (again, not shown for brevity) resemble the pre- 
vious example, the angular response is now much more 
selective (for both polarizations), and exhibits (in the 
P-polarization case) another peak for near-grazing inci- 
dence, which corresponds to a pseudo-Brewster angle of 
the ENG slab. The sensitivity to losses in the ENG ma- 
terial is illustrated in the inset of Fig. [H in terms of 
the peak transmittance as a function of the loss-tangent 
(at resonance), and qualitatively resembles the behavior 
observed in other tunneling phenomena involving SNG 
materials^ Finally, our simulations show that, for given 
parameters of the ENG slab, the sensitivity to variations 
within ~ 1/25 of the local wavelength with respect to the 
nominal DPS slab thicknesses ^2,3 would still yield mod- 
erate reductions (~ — 3dB) of the peak transmittance 
(with possible frequency shifts). 

To sum up, we have illustrated a counterintuitive EM- 
field resonant tunneling mechanism that can take place 
by pairing an SNG slab with a DPS bi-layer. In particu- 
lar, we have worked out analytically the design criteria, 
and explored the dependence on frequency, incidence di- 
rection, polarization, and material losses. Our results, 
which can be extended (along the lines of Refs. [HI and 
I20I ) to more general (e.g., optical, quantum-mechanical) 
asymmetrical tunnel barriers, also allow direct compar- 
isons with ENG-MNG paired configurations, and can be 
understood in terms of the equivalent wave-impedance 
properties exhibited by a DPS bi-layer and a "matched" 
(according to Ref. SNG slab. This observation may 
suggest more general application scenarios of such DPS 
bi-layer, for which SNG-like responses may be emulated 
via simpler dielectric slabs, which may somehow compen- 
sate, in a simple geometry, opaque ENG metamaterial 
slabs for a variety of applications in which complemen- 
tary metamaterials are paired together^ Also of inter- 
est is the possible adaptation for applications to super- 
resolution imaging schemes. 
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